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Effects of Refractive Index and 
Diffuse or Specular Boundaries on a 
Radiating Isothermal Layer 

R. Siegel 1,3 and C. M. Spuckler 2,3 


Nomenclature 

a - absorption coefficient of layer* m _l 
c = specific heat of solid, W • s/kg • K 
D = thickness of plane layer* m 
exb « blackbody spectral emissive power, 

W/m 2 • Atm 

£ 3 = exponential integral function* £3 (?) = Jo M exp 

F m = blackbody fraction of mth spectral band 
n * index of refraction 

q = heat flux, W/m 2 ; spectral flux, W/m • /tm 
= externally incident spectral radiation fluxes, 
W/m 2 • 

T e = temperature of environment, K; f, = 7>To 
7 0 = layer temperature at time / = 0, K 
T w = absolute temperature of layer, K; f w = T„/T 0 
T u T 2 = temperatures of radiating surroundings, K 
t - time, s; dimensionless time, t = toT^/pcD 
a, « = absorptance and emittance of plane layer 
B = angle of incidence 3 

p = reflectivity of interface; density of solid, kg/m 
a = Stefan-Boltzmann constant, W/m 2 • K 4 
r = transmittance of path at angle x> t = exp 
( - aD/ cos x) 

X = angle of refraction 

Subscripts 

dif = for diffuse surfaces 
/* o - incoming and outgoing radiation 
m - the mth spectral band 
5 = for specular surfaces 
X = spectral quantity 

j_ , R = components for perpendicular and parallel po- 
larization 

Superscripts 

/, 0 — at inside or outside surface of an interface 


Introduction 

Ceramics are under development for high-temperature ap- 
plications such as in combustion chambers of high-speed air- 
craft. Some ceramics are reinforced with ceramic fibers, and 
to protect the composite from the combustion environment it 
is coated with another ceramic to form a laminated layer. To 
model these materials thermally in order to study possible 
radiation effects arising from their partial transparency, it is 
necessary to specify the nature of reflections at their interfaces. 
A convenient assumption is that reflections are diffuse; this 
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Fig. 1 Plane layer with specular or diffuse boundaries and with diffuse 
or normally Incident radiation 


may be realistic because some external surfaces and internal 
interfaces of ceramics are rough. Some surfaces, however, may 
be smooth so they are specularly reflecting. This note will 
examine in a simple way by considering an isothermal layer, 
some effects of diffuse versus specular boundary conditions. 
In some cases the diffuse condition yields increased temper- 
atures due to augmented retention of incident energy by in- 
ternal reflections. This is a function of the refractive index of 
the material. Much of the work in the literature is for glass, 
n • 1.5; for ceramics n can be significantly larger, thus pro- 
viding increased internal radiation effects. The limiting con- 
dition is examined where there is no scattering inside the layer 
as this demonstrates the maximum effect of the surface con- 
dition. For specular surfaces, a small amount of scattering 
within the layer produces a partially diffuse behavior; for large 
scattering, diffuse behavior will result. 

Analysis 

Absorptance of Layer With Specular Boundaries and No 
Internal Scattering, The layer in Fig. 1 with refractive index 
n > 1 has incident diffuse or beam energy. The layer partially 
absorbs radiation, but does not scatter internally. For spec- 
ularly reflecting surfaces, radiation in each incident direction 
6 is refracted and internally reflected as given by Fresnel re- 
lations. The fraction of energy absorbed is given as (Siegel and 
Howell, 1992, p. 929): a(6) = [1 - P°(B) ] [1 - t( 0)]/[1 - 
p°(0)r(0)] where the internal transmittance is r(0) = exp [ - aD/ 
cos x(0)l and the angle of refraction is x(0) = sin 1 K sin W 
n]. The reflectivity p°(B) is a function of the polarization com- 
ponent so a(0) has values corresponding to the perpendicular 
and parallel components. For diffuse incidence of unpolarized 
radiation, the absorptance of a layer with specular surfaces is 
found by integrating the absorbed energy for all incident solid 


angles to give 

r/2 


<*,(«, aD)= j 

l<*j.(0) + ai(0)l cos B sin 6 dB 
0 

(1) 

where the p° needed to obtain a L and ai are 


<2> 

ii 

(n 2 -sin 2 0) 1/2 -cos 0~| 2 

(2a) 

(n 2 -sin 2 0) l/2 + cos 
n 2 cos B-(n 2 - sin 2 0) I/2 1 2 

P°M = 

(2b) 

n 2 cos 0 + (n 1 - sin 2 0) 1/2 J 
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With all quantities in Eq. (1) in terms of 0, the integration is 
performed numerically for various n and aD. For the special 
case of normal incidence, a s = (1 - p°) (1 - r)/(\ - p°r) 
where p° = [(/? - 1 )/(n + l)] 2 and r = exp(-aD). 

For emission from an isothermal layer, the local internal 
emission is considered for each direction. Radiation internally 
incident on a boundary is totally reflected if it arrives at an 
angle larger than the critical angle for reflection, which depends 
on n. The formulation by Gardon (1956) for the emittance of 
an isothermal layer with specular boundaries provides the same 
result as a, in Eq. (1); this is expected from Kirchhoff s law. 


Absorptance of a Layer With Diffuse Boundaries and No 
Internal Scattering. The layer boundaries are now considered 
diffuse. For a rough surface it is assumed that incident radia- 
tion will strike surface elements at random angles so the surface 
reflectivity is independent of the angular distribution of the 
incident energy. Hence all energy incident upon and leaving a 
boundary is considered to behave in a diffuse manner. The 
diffuse reflectivity is given by the hemispherically averaged 
relation of the Fresnel equations (Siegel and Howell, 1992, p. 
115), 


Pdir(n) 


1 (3n+ l)(/r- I) n 2 (n 2 - l) 2 

2 + 6(n+ l) 2 + (w 2 +l) 3 



2/i 3 (n z + 2/i- 1) 8w 4 (w 4 + 1) ■ 

(« 2 + 1)(« 4 - 1) + (n 2 + l)(/i 4 - l) 2 " * 


(3) 


The absorptance of the layer is obtained by considering its 
interaction with a unit incident flux, q f] = 1. The outgoing 
flux from the internal side of each diffuse boundary (see Fig. 
1) is then 

<7o( 0)« 1 -Pdif + <7/(0)p dif ; Q 0 (O) =qi(D)p l dlf (4a, b) 

where p dif is at the internal side of a diffuse boundary. As a 
result of transmission within the layer, the fluxes at the two 
internal sides of the boundaries are related by (Siegel and 
Howell, 1992, p. 707) 


qi(0) = 2q o (D)Ei(aD); q^D) =2q o (0)E>(aD) (5a, b) 


Equations (4) and (5) are solved for <? o (0), q 0 (D)> 9/(0), and 
q,(D ). The fraction of incident energy that is absorbed is a dif 
= 1 - (1 - piif) [ 9 /( 0 ) + qAD ) ], which yields 


Ofdir ( n » fl£>) - ( 1 “ Pdif) 


l-2£ 3 (aZ>) 

1 -pi if 2E 3 (flD) 


( 6 ) 


To place a dl f in terms of the reflectivity p d if on the outer surface 
of a boundary, the relation can be used that (Richmond, 1963) 
pi if = i - (I//* 2 ) (i - pit). 


Heal Balance Relations for Layer Temperature With Trans- 
fer Only by Radiation. The layer absorptances are spectral 
quantities since the absorption coefficient depends on wave- 
length. Two examples are given to demonstrate effects of dif- 
fuse and specular boundaries. The first is the equilibrium 
temperature attained by an isothermal layer subjected to ra- 
diation incident on both boundaries and with radiation being 
the only means of energy transfer. The balance of absorbed 
and emitted energy is 

1 (<?Ari + q\r 2 )ot\(n, a x D)dk 
J o 

= 2 f e x (n, a h D)e xb {\ t T w )d\ (7) 
J o 

where t x - cr x . Equation (7) is solved for the equilibrium 
temperature T w . 

A second example is the radiative cooling of a layer initially 
at temperature T l)t when placed in surroundings at T\ = T 2 
= 7*,. This is governed by the transient energy equation 


-pcD 


dT w , 

f- 2 ). 


axO)[ew,(X, T w ) 


-e xb (\, 7-,)]dX (8) 

For a spectrally banded calculation with M bands with each 
band designated by m, 

. dT w 


-pcD-^ = 2 2 amD )\F m (T v )aTi~ F m (T t )aT\\ 


(9) 


This is placed in dimensionless form and integrated to give 

•' dt : 


r ( 7V>=M s 

sr 


2 €m(«. a m D)[F m {f:T 0 )r:-F m a r )ft] 


( 10 ) 

Equation (10) is numerically integrated to yield cooling curves 
for various T 0 and T e . 


Results and Discussion 

Absorptance of Layer. For specular boundaries the layer 
absorptance was evaluated from Eq.(l) for both diffuse and 
normal incidence of external radiation. Results are in Fig. 2(a) 
for refractive indices from 1 to 4 and optical thicknesses up 
to 4. For a x £> larger than 4 almost all incident energy that is 
not reflected is absorbed in the layer. The layer absorptance 
then approaches an asymptotic value of one minus the external 
reflectivity of the boundary. Since the boundary reflectivity is 
smaller for normal incidence than for diffuse, the a Xt5 is larger 
for normal incidence. For a\D < 4 the a x > is influenced by 
the radiative path lengths within the layer. For a large n, 
incident diffuse energy is refracted into a small cone angle 
around the normal direction through the layer. Radiation within 
the nonscattering layer then travels along paths not much longer 
than the path length D for normal incidence; hence, for large 
n the a X4 becomes independent of normal or diffuse incidence. 
There are differences in a Xt5 when n is near 1 and refraction 
is therefore small. Diffuse radiation then travels along longer 
paths and a XtJ becomes increased. This is compensated some- 
what by the reflectivity of the external boundary being larger 
for diffuse incident radiation than for radiation incident in the 
normal direction. 

In Fig. 2(b) the absorptance of a layer with diffuse bound- 
aries is compared with that for specular boundaries; the in- 
cident radiation is diffuse. The nature of the boundary has a 
very significant effect when Q\D is less than about 3. For n = 

1 there are no surface reflections so the results are the same 
for the solid and dashed curves. As n is increased the absorp- 
tance values become increasingly different, so for n = 4 with 
a x D < 3 there is a substantial difference in ot\. For specular 
boundaries a x decreases as the angle of refraction becomes 
smaller with larger n t thus providing shorter path lengths for 
absorption in the layer. For diffuse boundaries the increase in 
a x with n is attributed to radiation within the layer being 
diffused at the boundaries into directions larger than the critical 
angle for reflection. This tends to retain energy within the layer 
and increase its absorption as radiation travels internally by 
multiple reflections. The effect that this can have on the layer 
temperature is shown by results that follow as calculated from 
Eqs. (7) and (10). If a layer with specular boundaries has 
internal scattering it will act more like a diffuse layer; the 
results here are without internal scattering and are expected to 
show the maximum possible effect of the two different types 
of boundaries. 

Effect of Surface Condition on Layer Temperature. To 

illustrate how surface condition can affect the layer temper- 
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TEMPERATURE OF HOT RESERVOIR, T, , K 



REFRACTIVE INDEX, n 

Fig. 2 The effect of diffuse or specular boundaries, and beam or diffuse 
Incidence, on spectral absorptance (or spectral emittance) of a plane 
layer without Internal scattering as a function of optical thickness and 
refractive index: (a) Boundaries are specular and incident external ra- 
diation Is diffuse or Is In the normal direction; ( b ) boundaries are diffuse, 
or are specular with diffuse incident radiation 


ature, results are given for a layer with D = 0.5 cm and a x 
approximated as constant in each of three spectral bands, 
similar to the properties of glass: a x = 0.2, 4, 100 cm' 1 for 
A = 0-2.7, 2. 7-4. 4, 4.4-oo ^m. Using these properties the 
spectral absorptance and emittance in each band was obtained 
from Eqs. (1) and (6). The equilibrium temperature T w for 
heat transfer only by radiation was then obtained from Eq. 
(7), and transient radiative cooling results were evaluated from 
Eq.(10). 

Figures 3(a) and 3(6) show how equilibrium temperatures 
are affected for two different conditions of the incident ra- 
diation. As shown by Fig. 2(6) , when a x D is small the specular 
and diffuse boundary conditions produce significantly differ- 
ent a x . Hence for a layer without internal scattering the con- 
dition of the boundary can produce large differences in energy 
absorption if there is large incident energy in wavelength re- 
gions where a x D is small. For a high temperature radiant source 



BEAM HEAT FLUX, kW/m 2 

Fig. 3 Effect of diffuse or specular boundaries on the equilibrium tem- 
perature of a layer with a refractive index of 1.5 or 3, and for the limiting 
case n = 1. The spectral absorption properties are similar to glass (a x 
= 0.2, 4, 100 cm for X » 0-2.7, 2.7-4.4, 4.4-ao j*m) and there is no 
internal scattering; T 3 = 300 K. (a) Incident radiation on one side Is from 
a hot reservoir at T|. (b) Incident radiation on one side is a beam normal 
to the layer with a spectral distribution like that of solar radiation. 


the incident energy fraction is large in the short wavelength 
region where semitransparent materials can have a small a x . 
The wavelength selective heat balance can then be sensitive to 
the specification of specular or diffuse boundary conditions. 
Reradiation is at longer wavelengths where a x D is generally 
large; hence from the results in Fig. 2(6) the reradiated energy 
is not affected significantly by the boundary being specular or 
diffuse. For a moderate source temperature both the incident 
and emitted energy will be at long wavelengths where a x D can 
be large causing a x to be insensitive to specular or diffuse 
boundaries. 

In Fig. 3 (a) the layer temperature T w is given for radiative 
heating by blackbody reservoirs with T { given on the abscissa, 
and 7\ = 300 K. The T w is given for diffuse and specular 
boundaries with the layer refractive index /? = 1, 1.5, and 3. 
For n = 1 there are no surface reflections so the energy balance 
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DIMENSIONLESS TIME, taTj/pcD 


Fig. 4 Transient cooling of a layer with Initial temperature 7 0 = 1000 
K (thin lines) or 2000 K (thick lines), T # = 300 K, and with the same 
absorption properties as for Fig. 3 (a x = 0.2, 4, 100 cm * 1 for \ * 0-2.7, 
2.7-4.4, 4.4-oo fi m); comparisons are made for diffuse and specular 
boundaries with n = 1,5 and 3 (for n * 1 there are no surface reflections) 

depends only on internal absorption and emission. For n > 1 
the higher a\ for small a\D with diffuse boundaries produces 
increased absorption. The resulting T w for diffuse boundaries 
is somewhat higher than for the specular case and the difference 
increases to a maximum of about 50 K as n and T\ are increased. 
The temperature for an opaque black layer (or an opaque gray 
layer with diffuse boundaries) is shown as an upper limit. From 
Eq. (7) this is given by T w = [(T, 4 + r 2 4 )/2] l/4 . 

For Fig. 3(b) the nature of the incident radiation is changed 
on the hot side of the layer. The incident energy is a beam 
normal to the layer with its heat flux given on the abscissa. 
Two limiting cases are given for n — 1 ; the solid curve is for 
the beam being diffused at the boundary, and the dot-dash 
curve is for the beam remaining normal to the boundary as it 
passes into the layer. The heat flux is high (for comparison 
the solar constant at earth orbit is about 1 .35 kW/m 2 ) and has 
the same spectral distribution as solar energy (black body spec- 
trum at 5780 K). This spectrum provides a much larger energy 
fraction in the short wavelength region than the hot reservoir 
in Fig. 3(a) and produces substantially higher T w when the 
boundaries are diffuse. The upper limit is shown for a black 
opaque plate as given by T w = + Ti)/2 ] WA . 

Effect of Surface Condition on Transient Radiative Cool- 
ing. Transient radiative cooling results from Eq. (10) are in 
Fig. 4 for two initial temperatures T 0 = 1000 K (thin lines) 
and 2000 K (thick lines), and T e = 300 K. The absorption 
properties are the same as for Fig. 3. For n = 1 there are no 
surface reflections and hence no effect on cooling of the bound- 
ary being specular or diffuse. For T 0 = 1000 K most of the 
energy radiated by the layer is in the long wavelength region 
where a\D is large for the spectral properties used here. For 
large a^D, Fig. 2(b) shows that increasing n decreases the ax 
= * x so the cooling is less rapid. The is smaller for specular 
surfaces, which produce longer cooling times. The situation 
changes somewhat when T 0 is increased to 2000 K. A greater 
fraction of the energy is then in the short wavelength regions 
where a x D is small. This changes the trend for the solid curves 
(diffuse surfaces) in Fig. 2(b) which, for small a x D, show ex 
increasing with n to a maximum. This enhances the cooling 
rate with the increase in n as shown in Fig. 4 by the thick solid 
curves as compared with the n = 1 results (thick dashed curve). 


The limiting results are also shown for cooling of a black 
opaque layer. In this instance there is no effect of spectral 
selectivity and the layer with highest initial temperature T 0 
cools most rapidly. 

Conclusions 

Equilibrium temperatures of an absorbing-emitting layer 
were obtained for exposure to incident radiation and with the 
layer boundaries either specular or diffuse. For high refractive 
indices the surface condition can influence the radiative heat 
balance if the layer optical thickness is small. Hence for a 
spectrally varying absorption coefficient the layer temperature 
is affected if there is significant radiative energy in the spectral 
range with a small absorption coefficient. Similar behavior 
was obtained for transient radiative cooling of a layer where 
the results are affected by the initial temperature and hence 
the fraction of energy radiated in the short wavelength region 
where the absorption coefficient is small. The results are for 
a layer without internal scattering. If internal scattering is 
significant, the radiation reaching the internal surface of a 
boundary is diffused and the effect of the two different surface 
conditions would become small. 
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Nomenclature 

h - surface convection heat transfer coefficient 
H - channel width = r H - r 0 
K - thermal conductivity 
m - fluid mass flow rate 
Nu = local Nusselt number = hr 0 /K y 
r, x = coordinate distances in the radial and axial 
directions, respectively 
r 0 = radius of the moving fiber 
R, X = dimensionless r and x, R = x/r 0 , X = x/r 0 
T - local physical temperature 
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